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' Abstract. Revised upper limits on gamma-ray burst repetition rates are found using 

(N ■ the BATSE 3B and 4B catalogs. A statistical repetition model is assumed in which 

sources burst at a mean rate but in which BATSE observes bursts randomly from each 
source. 



O : INTRODUCTION 
+^ 

Regardless of the gamma-ray burst distance scale, the question of whether 
bursters repeat (e.g. [8] [11]) is of great importance to our understanding of 
^ I them: their energy production mechanisms are constrained by repetition. There 

c3 ' is strong statistical evidence to indicate that bursters are not observed to repeat 

on timescales of days to months ( [9] [4] [2]), but the recent identification of four 
gamma-ray bursts close together in both location and time [3] has left this possi- 
bility open for at least some bursters. Of equal interest is the possibility that they 
might repeat on timescales of years or longer, and that repetition might take some 
form where they do not burst steadily. 

Detection of repetition is complicated by burst localization errors (do two bursts 
belong to the same source?) and by incomplete angular sky exposure and hmited 
trigger efficiency (what fraction of source repetitions is observed?). The repeti- 
tion question can be studied using two observational parameters [6] [4]; / and (z/). 
The variable / represents the fraction of sources that produce more than one de- 
tected burst, while (z^) represents the mean number of detected bursts per source 
producing more than one detected burst. 



ANALYSIS 



Localization Errors. Localization errors are modeled as a combination in 
quadrature of a 1.6° systematic error and a statistical error. The statistical error 
is a combination of a fiuence-dependent mean and a fluence-dependent random 
component obtained from the BATSE 4B Catalog. 

Sky Exposure. The sky exposure used is that hsted in the BATSE 4B Catalog. 
It should be noted that this exposure is only valid for non-overwriting bursts, so 
that overwriting bursts must be removed from the dataset. 

Trigger Efficiency. The effects of trigger efficiency on repetition are not in- 
cluded in this analysis. This is due in part to calibration work still in progress 
(see Pendleton, Hakkila, & Meegan, this conference) and in part to the model- 
dependent way in which the repetition luminosity function enters into the analysis 
(e.g. [1] [4]). 

The Burst Catalogs. For this analysis we choose to independently analyze 
the BATSE 3B [7] and 4B catalogs. The catalogs have had all overwriting bursts 
removed, leaving 1060 3B bursts and 1554 4B bursts. We recognize that 4B bursts 
have been detected using several triggering energy channel combinations; this com- 
plicates the analysis in complex ways not modeled here. The 3B analysis does not 
suffer from this effect. 

Clustering Statistic. The value of the Two-Point Angular Correlation Func- 
tion found within 7.2° [TPACF(6' < 7.2°)] is used as the clustering test for this 
analysis. This angle has been found to contain an optimum signal for repeating 
bursts [10] due to the overall distribution of burst localization probabilities. 

Repetition Models. Two models arc used; a Euclidean model and a model 
based on an Einstein-de Sitter cosmology with the faintest bursts at redshift z — 2 
(Hereafter referred to as the ES2 model). Both spatial models have only one free 
parameter: it is assumed that bursters burst with a mean intrinsic repetition rate 
-Rint (repetitions source^^ ycar^^), but that equal probabilities exist that any burster 
might be the next one to burst. By varying i?int, various fractions of the bursts 
/ are found to belong to clusters composed of various numbers of detected bursts 
V (having an average value (z^)), which corresponds to a mean observed repetition 
rate i?obs. In the limit of very low mean intrinsic repetition rates (i?int ~ 0), these 
models are essentially the same as non-repeating models, since no clusters of two 
or more bursts will be produced during the BATSE operating time. In Euclidean 
space, i?int is always less than i?obs; since the overall exposure is less than unity. In 
the case of the ES2 model, i?int is assumed to be constant in the comoving frame. 
More distant ES2 bursters appear to repeat at a slower rate due to time dilation, 
so that a higher i?int is needed to produce a Euclidean equivalent i?obs. 

Monte Ceu"Io Simulations. Monte Carlo simulations are performed while vary- 
ing the free parameter i?int. It is found from the Monte Carlo simulations that 



-Robs — 0.481i?int ; Euclidean 

— 0.310i?int;ES2 



(1) 



Thus, there is no real difference between the Euchdean and ES2 cosmological 
models in terms of the repetition observables. 

The Monte Carlo simulations allow (v) vs / plots to be produced for both the 
BATSE 3B and 4B catalogs. These are presented in Figure 1, along with the corre- 
sponding mean observed repetition rates. i?obs from (u) and / differs between the 
two catalogs because the exposure differs. The 4B exposure is slightly better than 
the 3B exposure; if the two were identical the rates would coincide on the plots. If 
the rates were identical, then the longer 4B Catalog would produce more chistcr 
composed of larger numbers of bursts. This effect is enhanced by the better 4B 
exposure, so that even more and larger clusters should be produced. It is therefore 
expected that repetition rates of the type described should be easier to identify in 
the 4B Catalog than in the 3B Catalog. 

Repetition Limit Procedure. The BATSE 3B and 4B values of the TPACF 
within 7.2° of one another [TPACF(^ < 7.2°)] are slightly negative (the 3B value 
of —0.070 is 1.5(7 from the isotropic value of zero, while the 4B Catalog value of 
—0.035 is 1.1(7 from isotropy), indicating a depletion of close burst pairs. Since 
models with nonzero repetition rates produce positive values of this statistic, both 
BATSE catalogs are most consistent with no repetition. 

We can also ask how consistent the results are with different mean observed rep- 
etition rates i?obs, so that we can get some idea concerning upper limits allowed on 
repetition rates. To do this, the TPACF(6' < 7.2°) measurement (with associated 
error distribution; assumed Gaussian) at each i?obs can be compared to the mini- 
mum expected value of TPACF < 7.2°) = at i?obs = (assuming that TPACF 
= is the correct value). This sensitivity calculation can provide some idea as to 
the maximum rate i?obs that would go undetected by the BATSE instrument. 
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FIGURE 1. Correspondence between the mean observed repetition rate i?obs and the Meegan 
/ and {v) parameters for the repetition modehn question. Specific rates Rohs (bursts source"^ 
year"-^) are indicated for both the BATSE 3B (open circles; itahc lettering) and 4B (closed circles; 
roman lettering) catalogs. The shape of the {v) vs. / relationships seen here is characteristic of 
this repetition model. 



TABLE 1. Limits on mean observed and intrinsic repetition 
rates (in units of bursts source"^ year~^) for Euclidean and 
ES2 models. Limits indicate rates for which the probability is 
90% that TPACF(6' < 7.2°) > 0, as opposed to the negative 
TPACF(6I < 7.2°) values obtained from the BATSE 3B and 4B 
catalog.s. 



Catalog 


i?obs Limit 


-Rint; Euclidean Limit 


•Rint;RS2 Limit 


BATSE 3B 


0.15 


0.31 


0.48 


BATSE 4B 


0.05 


0.10 


0.16 



CONCLUSIONS 

Figure 2 indicates the probability that TPACF(^ < 7.2°) > assuming different 
i?obs rates in the BATSE 3B and BATSE 4B catalogs. At the 90% confidence level, 
-Robs < 0.15 bursts source" year~^. The larger data set of the BATSE 4B Catalog 
would place even greater constraints on the observed repetition rate, if one considers 
the different triggering criteria used in this catalog to be unimportant. For the 4B 
Catalog, it is unlikely that i?obs > 0.04 bursts source"^ year~^. This indicates that 
repetition is either absent or that sources only repeat rarely on the timescales of the 
BATSE catalogs. These 90% probabilities correspond to model-dependent limits 
on the mean intrinsic repetition rates, which are indicated in Table 1. 

In a previous work, we [4] quoted similar types of limits on the Euclidean mean 
intrinsic repetition rate, finding that i?obs > 0.05 bursts source"^ year~^ was un- 
likely in the BATSE 3B Catalog. Our revised 3B 90% upper limit is higher at 
^obs — 0.31 bursts source"^ year~^, in part due to a more direct definition of our 
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FIGURE 2. Probability of obtaining a value of TPACF(6' < 7.2°) > vs. mean observed 
repetition rate i?obs- The plot indicates that the larger 4B data set (closed circles) is more sensitive 
to detecting repetition of the type described, than the 3B data set (open circles). Unknown effects 
in the 4B, such as that of triggering on different energy channels, have not been taken into account 
in this analysis. 
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statistical significance, and in part indicating that BATSE is not as sensitive to 
detecting repetition as a result of the improved sky exposure analysis. The new- 
sky exposure requires that overwrites be excluded; with fewer bursts, higher rep- 
etition rates are needed to be detected above the random background. The new 
exposure is also higher than that estimated before, indicating that more of any 
intrinsic repetitions should be detected. By selecting fewer bursts from sources, 
cluster sizes tend to be smaller. Thus [i/) is smaller for a given repetition rate, and 
is more difficult to detect. 

Our analysis technique could still be more sensitive to detecting repetition than a 
recent technique that combines localization probabilities from multiple experiments 
to obtain smaller localization regions [5]. That technique limits its database to 
bursts detected by both experiments, so that the sky exposure and trigger efficiency 
for the combined experiment must then be found by multiplying these detection 
probabilities together. We have shown here that limits on repetition rates depend 
both on localization uncertainty and sky exposure. 

Our repetition model allows that some observed clusters can be large, and in fact 
suggests that is it possible for the October 1996 burst "cluster" [3] to have come 
from one source. From the upper limit on the repetition rate described above, it is 
expected that at most 29 clusters of 3 bursts and 4 clusters of 4 bursts should be 
observed in the BATSE 3B Catalog. The tighter hmits on 4B Catalog repetition, 
however, indicate that only as many as 19 clusters of 3 bursts and 1 cluster of 4 
bursts are expected. 
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